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ABSTRACT 


A  comprehensive  data  set  from  the  ocean  and  atmosphere  was  obtained  just  north 
of  the  Monterey  Bay  as  part  of  the  Monterey  Bay  2006  (MB06)  field  experiment.  The 
wind  stress,  heat  fluxes,  and  sea  surface  temperature  were  sampled  by  the  Naval 
Postgraduate  School’s  TWIN  OTTER  research  aircraft.  In  situ  data  were  collected  using 
ships,  moorings,  gliders  and  AUVs.  Four  data-assimilating  numerical  models  were 
additionally  run,  including  the  Coupled  Ocean/ Atmosphere  Mesoscale  Prediction  System 
(COAMPS®)  model  for  the  atmosphere  and  the  Harvard  Ocean  Prediction  System 
(HOPS),  the  Regional  Ocean  Modeling  System  (ROMS),  and  the  Navy  Coastal  Ocean 
Model  (NCOM)  for  the  ocean. 

The  scientific  focus  of  the  Adaptive  Sampling  and  Prediction  Experiment  (ASAP) 
was  on  the  upwelling  /  relaxation  cycle  and  the  resulting  three-dimensional  coastal 
circulation  near  a  coastal  promontory,  in  this  case  Point  Aho  Nuevo,  CA.  The  emphasis 
of  this  study  is  on  the  circulation  over  the  continental  shelf  as  estimated  from  the  wind 
forcing,  two  ADCP  moorings,  and  model  outputs.  The  wind  stress  during  August  2006 
consisted  of  3-10  day  upwelling  favorable  events  separated  by  brief  1-3  day  relaxations. 
During  the  first  two  weeks  there  was  some  correlation  between  local  winds  and  currents 
and  the  three  models’  capability  to  reproduce  the  events.  During  the  last  two  weeks, 
largely  equatorward  surface  wind  stress  forced  the  sea  surface  and  barotropic  poleward 
flow  occurred  over  the  shelf,  reducing  model  skill  at  predicting  the  circulation.  The 
poleward  flow  was  apparently  remotely-forced  by  mesoscale  eddies  and  alongshore 
pressure  gradients,  which  were  not  well  simulated  by  the  models.  The  small,  high- 
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resolution  model  domains  were  highly  reliant  on  correct  open  boundary  conditions  to 
drive  these  larger-scale  poleward  flows.  Multiply-nested  models  were  no  more  effective 
than  well-initialized  local  models  in  this  respect. 
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1.  Introduction 


The  Adaptive  Sampling  and  Prediction  Experiment  (ASAP)  was  launched  during 
the  summer  of  2006  with  the  over  all  goal  to  learn  how  to  deploy,  direct,  and  utilize 
autonomous  vehicles  most  efficiently  to  sample  the  ocean,  assimilate  the  data  into 
numerical  models  in  real  or  near-real  time,  and  predict  future  conditions  with  minimal 
error.  A  fleet  of  four  Spray  and  nine  Slocum  gliders  [Davis  et  ah,  2002]  under 
autonomous  control  were  deployed  in  late  July  to  cooperatively  sample  a  control  volume 
in  an  upwelling  center  for  a  period  of  30  days.  The  scientific  focus  was  to  close  the  heat 
budget  in  a  three-dimensional  upwelling  center  and  thus  illuminate  the  along-  and  across- 
shore  transport  and  fundamental  dynamics  of  the  process.  The  glider  fleet  patrolled  the 
boundaries  of  the  control  volume  and  a  research  aircraft  observed  the  fluxes  through  the 
sea  surface.  Two  bottom-mounted  acoustic  Doppler  current  profilers  (ADCPs)  were 
deployed  on  the  54  and  92  m  isobaths  to  monitor  the  currents  over  the  continental  shelf 
and  provide  ground  truth  for  the  newer,  more  innovative  instrumentation. 

Based  on  previous  work,  the  evolution  of  coastal  upwelling  near  Point  Ano 
Nuevo  can  be  described  as  follows:  By  comparing  aircraft-sensed  winds  and  SST,  Ramp 
et  al.  [2009]  detennined  that  cold  water  <  13°C  first  appears  at  the  Point  within  12-24 
hours  of  the  onset  of  upwelling-favorable  winds.  As  the  wind  continues  to  blow  for  5-10 
days,  a  cold  plume  forms  which  may  move  further  offshore  or  turn  southward  across  the 
Monterey  Submarine  Canyon  (MSC)  [Rosenfeld  et  al.,  1994,  Ramp  et  al.,  2005;  Haley  et 
al.,  2009].  A  warm,  fresh,  anticyclonic  feature  (the  “Monterey  Bay  Eddy”  or  MBE)  often 
exists  offshore,  which  interacts  with  the  cold  plume  in  some  complex  way  [Rosenfeld  et 
al.,  1994;  Paduan  and  Rosenfeld,  1996;  Ramp  et  al.,  2005].  It  has  now  been  determined 
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during  the  second  Autonomous  Ocean  Sensing  Network  (AOSN-II)  experiment  that  the 
position  and  strength  of  the  cyclonic  feature  can  detennine  the  bifurcation  of  the 
upwelling  front  [Ramp  et  ah,  2009].  When  the  MBE  lies  offshore,  the  upwelling  plume 
is  transported  offshore  by  Ekman  advection.  When  the  feature  has  moved  inshore,  the 
combination  of  blocking  and  geostrophic  current  force  the  plume  to  move  southward 
across  the  mouth  of  the  bay.  Time  series  of  sea  surface  temperature  from  aircraft  and 
AVHRR  images  show  that  the  MBE  moves  rapidly  onshore  during  wind  relaxation 
events  and  retreats  back  offshore  when  the  winds  re-intensify  [Rosenfeld  et  ah,  1994; 
Ramp  et  al.,  2005].  Moored  observations  however  have  also  shown  an  onshore  and 
southward  translation  of  the  MBE  without  a  wind  relaxation  [Ramp  et  al.,  1997].  The 
mechanism  for  this  across-shore  eddy  translation  is  thus  not  fully  understood.  How  far 
the  eddy  forcing  penetrates  up  onto  the  continental  shelf  has  also  not  been  well  observed. 
Since  these  features  are  more  than  1000  m  deep  [Ramp  et  al.,  1997]  one  expects  that  due 
to  rotational  rigidity,  their  penetration  onto  the  shelf  would  be  limited,  similar  to  Gulf 
Stream  Warm  Core  Rings  [Ramp  et  al.,  1989]. 

Much  less  obvious  is  the  ocean  dynamics  during  the  wind  relaxation  events, 
defined  here  to  be  periods  of  poleward  or  no  wind.  During  these  events,  barotropic 
poleward  flows  typically  appear  over  the  continental  shelf  along  the  west  coast  of  the 
U.S.  Observational  examples  of  this  include  the  Coastal  Ocean  Dynamics  Experiment 
(CODE)  region  off  Point  Arena,  CA  [Winant  et  al.,  1987;  Lentz,  1987],  further  south  off 
Point  Sur,  CA  [Chelton  1984;  Ramp  and  Abbott,  1998],  and  south  of  Cape  Blanco,  OR 
[Ramp  et  al.,  2008],  Recent  modeling  studies  attribute  this  current  to  a  poleward 
alongshore  pressure  gradient  force  (PGF)  south  of  the  major  headlands,  which 
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encourages  both  poleward  flow  near  the  coast  and  onshore  transport  in  geostrophic 
balance  with  the  poleward  PGF  [Gan  and  Allen,  2002a,  2002b].  Since  these  gradients  are 
typically  remotely  forced  over  very  large  scales,  there  is  usually  less  coherence  between 
local  winds  and  currents  during  relaxation  events  than  during  upwelling  events,  making 
them  much  more  difficult  to  reproduce  with  local,  high-resolution  numerical  models.  In 
fact,  uncertainties  of  coherent  structures  estimated  by  model  simulations  (Lermusiaux, 
2006;  Lennusiaux  et  ah,  2006)  during  such  relaxation  periods  have  large  uncertainties. 

In  this  paper,  a  rigorous  description  of  the  currents  over  the  continental  shelf  just 
south  of  Point  Ano  Nuevo  and  their  response  to  local  and  remote  forcing  is  presented 
using  in-situ  data  from  moored  ADCPs,  remotely-sensed  data  from  a  low-flying  aircraft, 
and  model  output  from  three  numerical  models  run  during  the  ASAP  experiment.  The 
paper  is  intended  to  set  the  stage  for  a  more  detailed  dynamical  analysis  of  the  three- 
dimensional  circulation  in  the  ASAP  region  using  all  available  assets.  Following  a  basic 
description  of  the  observational  tools  and  techniques  in  Section  2,  the  observed  winds  and 
currents  are  described  in  Section  3  and  compared  using  correlation  and  spectral  analysis. 
The  three  models,  including  their  grids,  initialization  and  data  assimilation  techniques, 
and  computational  schemes  are  set  forth  in  Section  4,  where  the  model  outputs  are  also 
quantitatively  compared  with  the  observed  currents  at  the  two  ADCP  locations.  The 
results  are  discussed  in  Section  5  and  summarized  in  Section  6. 


2.  Data  and  Methods 
2,1  Real-Time  ADCP  Data 
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The  moorings  deployed  for  the  August  2006  experiment  consisted  of  two  300  kHz 
acoustic  Doppler  current  profilers  (ADCPs)  in  trawl-resistant  bottom  mounts  (TRBMs) 
located  6.55  km  apart.  ADCP  1  was  located  at  36°  55.336'N,  122°  07.344'W,  54  m 
depth,  and  ADCP  2  was  located  at  36°  53.683'N,  122°  1 1.244'W,  92  m  depth  (Figure  1). 
The  data  were  sampled  in  4-m  bins  from  the  bottom  to  the  surface,  however  the 
uppermost  3  bins  (12  m)  were  lost  due  to  side-lobe  reflection.  Bottom-mounted  moorings 
are  advantageous  because  they  have  no  extraneous  movement,  and  allow  faster  sampling 
and  less  averaging  to  achieve  a  good  signal-to-noise  ratio  than  a  surface  mooring.  The 
basic  1 -minute  time  series  were  obtained  by  ensemble  averaging  60  1-s  pings  in  the 
instrument.  Subsequent  averaging  and  filtering  to  separate  the  frequency  bands  of 
interest  were  accomplished  in  the  laboratory.  An  inverse-Fourier  truncation  method  was 
used  to  remove  internal  waves  with  periods  of  five  hours  or  less.  The  residual  time  series 
were  filtered  again  with  the  half-power  point  at  33  hours  [Beardsley  et  ah,  1985]  to 
separate  the  tidal  and  inertial  motions  from  the  wind  forced,  mesocale,  and  mean 
currents.  A  principal  axis  technique  was  used  to  rotate  the  u  and  v  components  into 
across-  and  alongshelf  components.  This  method  of  rotation  ensures  that  a  minimum  of 
the  alongshelf  flow  will  be  translated  into  the  across-shelf  flow  [Kosro,  1987].  The 
principal  axis  rotation  agreed  with  the  local  orientation  of  the  bathymetry  and  yielded 
angles  of  -58°  for  ADCP  1  and  -48°  for  ADCP  2,  with  respect  to  true  north. 

A  novel  aspect  of  the  experiment  was  that  the  data  were  transmitted  back  to  the 
laboratory  in  real-time  using  Benthos  Telesonar  acoustic  modems  operating  on  the 
Seaweb  undersea  network  [reference].  This  system  allowed  the  ADCPs  to  be  bottom- 
mounted  in  TRBMs,  completely  free  of  surface  motion,  but  still  able  to  transmit  the  data 
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to  the  laboratory  in  real  time.  The  ADCPs  were  programmed  to  send  back  the  most 
recent  estimate  via  satellite  link,  every  10  minutes.  A  redundant  system  with  two  surface 
gateway  buoys  and  three  intermediate  relay  points  was  used  to  ensure  reliability  and  keep 
the  acoustic  power  output  levels  low.  The  data  were  moved  from  the  gateway  buoys  to 
the  Naval  Postgraduate  School  via  the  iridium  satellite  communications  network.  The 
system  worked  flawlessly  for  18  days  but  then  failed  due  to  a)  sea  water  incursions  into 
the  iridium  antennas  and  b)  low  power  on  the  surface  buoys.  The  real-time  current  data 
were  used  to  understand  and  adapt  the  sampling  paths  of  the  ASAP  mobile  assets. 


2.2  Wind  Observations 

The  wind  time  series  was  recorded  by  buoy  M2  located  at  36°  42’  N,  122°  23’  W 
at  the  mouth  of  the  Monterey  Bay  and  maintained  by  the  Monterey  Bay  Aquarium 
Research  Institute  (MBARI).  Buoy  M2  was  slighter  farther  away  from  the  two  ADCPs 
than  similar  buoy  Ml,  but  was  less  sheltered  by  coastal  topography  and  therefore  more 
representative  of  the  offshore  winds.  The  wind  time  series  were  low-pass  filtered  with 
the  same  3  3 -hour  filter  as  the  ADCP  data  and  rotated  22°  into  the  along-  and  across-shore 
components,  aligned  with  the  coastal  mountain  range.  This  is  a  different  rotation  than 
the  currents  were  rotated  to  match  the  local  bottom  topography. 


2.3  Aircraft  Data 

The  airborne  measurements  during  ASAP  were  obtained  using  the  TWIN  OTTER 
aircraft  owned  and  operated  by  the  Center  for  Interdisciplinary  Remotely-Piloted  Aircraft 
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Studies  (CIRPAS)  at  the  Naval  Postgraduate  School.  Thirteen  flights  were  conducted 
daily  between  August  1  and  15,  2006  (except  August  6  and  13)  at  an  altitude  of 
nominally  33  m  above  the  sea  surface.  The  low  altitude  was  chosen  to  pass  beneath  the 
stratus  cloud  deck  that  often  covers  the  Monterey  Bay  in  summer,  and  to  minimize  the 
vertical  divergence  in  the  heat  and  momentum  flux  observations.  The  flight  path  covered 
the  area  from  the  center  of  the  bay  north  to  about  Pigeon  Point,  out  to  35  km  offshore. 
This  took  about  2.5  h  to  execute  beginning  usually  about  10  am  local  time,  subject  to 
weather  conditions.  To  collect  air/sea  fluxes  at  high  temporal  and  spatial  resolution,  the 
sample  area  focused  on  the  region  surrounding  Point  Ano  Nuevo  where  the  upwelling 
plume  was  first  developing.  The  atmospheric  parameters  sampled  by  the  aircraft 
included  air  temperature  and  dew  point,  atmospheric  pressure,  wind  speed  and  direction, 
the  turbulent  fluxes  of  heat  and  momentum,  and  total  particle  number,  as  well  as  aerosol 
and  cloud/fog  droplet  size  distributions.  In  the  ocean,  sea  surface  temperature  (SST), 
ocean-leaving  radiance  at  193  wavelengths,  sun  glint,  and  surface  roughness  were 
observed.  The  SST,  air  temperature  and  dew  point,  and  wind  data  were  made  available 
for  model  assimilation  in  near-real  time,  while  the  other  parameters  required  additional 
post-processing. 


3.  Observational  Results 
3.1  Atmospheric  Forcing 

August  2006  was  a  fairly  typical  summer  month  off  central  California  with 
periods  of  sustained  upwelling  favorable  winds  at  about  10  m  s'1  separated  by  brief 
periods  of  weak  poleward  or  no  winds  (Figure  2).  At  buoy  M2,  upwelling-favorable 
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wind  events  with  approximately  equal  strength  took  place  during  August  1-3,  8-12,  16- 
25,  and  28-31.  Poleward  winds  less  than  2.5  m  s'1  blew  during  August  4-7,  13-15,  and 
26-27.  A  unique  strength  of  the  ASAP  data  set  was  the  spatial  distribution  of  the  wind 
stress  as  sampled  by  the  TWIN  OTTER  aircraft.  The  wind  stress  was  computed  from  the 
bulk  aerodynamic  flux  fonnulae  using  the  wind  speed,  air  temperature,  relative  humidity, 
air  pressure,  sea  surface  temperature,  and  altitude  as  inputs  [Fairall  et  ah,  1996].  A 
turbulent  boundary  layer  was  assumed  and  Monin-Obukov  similarity  theory  was  used  to 
adjust  the  stress  down  to  10  m  height.  The  assumption  was  validated  by  the  bulk 
Richardson  Number  (R,)  which  was  small  and/or  negative  throughout,  indicating  a  fully 
turbulent  boundary  layer.  The  only  exception  occurred  on  August  5,  an  exceptionally 
calm  day,  but  R,  was  still  less  than  10  and  the  same  method  was  used  to  compute  the 
stress. 

A  selection  of  8  of  the  15  flights,  4  during  a  relaxation  event  and  four  during 
upwelling,  show  significant  sub-mesoscale  variability  in  the  offshore  wind  stress 
distributions  (Figure  3).  The  first  wind  relaxation  event  actually  lasted  longer  near  Point 
Ano  Nuevo  than  indicated  by  the  time  series  at  Buoy  M2,  and  spanned  6  days  from  2-7 
August.  Especially  noteworthy  during  this  time  was  the  wind  stress  maximum  located 
just  to  the  north  of  the  Point.  This  feature  looks  like  an  atmospheric  expansion  fan 
similar  to  those  previously  observed  in  the  lee  (to  the  south)  of  major  capes  during 
upwelling-favorable  winds  [Dorman  et  ah,  1995,  2000;  Gan  and  Allen,  2004;  Gan  et  ah, 
2004],  however  such  features  have  not  been  directly  observed  during  poleward  wind 
events.  The  maximum  wind  stresses  in  these  features  (0.08  N  m"2)  on  August  3ld  and  4th 
were  actually  much  stronger  than  indicated  by  the  offshore  buoys  and  were  counter  to  the 


11 


generally  held  idea  that  the  winds  are  strongest  offshore  and  weaken  as  the  coast  is 
approached.  The  weakest  winds  during  the  relaxations  were  near  the  mouth  of  the 
Monterey  Bay. 

During  the  upwelling  event  starting  on  August  8th,  all  aircraft  observations 
showed  an  alongshelf  wind  stress  maximum  off  Point  Ano  Nuevo  (Figure  3).  This 
maximum  was  strongest  on  August  8th  and  1 1th  (>0.16  N  in'2)  and  slightly  weaker  on  the 
9th  and  10th.  This  clear  maximum  in  the  alongshore  wind  stress  may  contribute  to  the 
formation  of  the  upwelling  center  off  Ano  Nuevo.  Consistent  with  previous  studies,  there 
was  a  wind  shadow  behind  the  coastal  mountains  in  the  northeast  comer  of  the  Monterey 
Bay  (Graham  and  Largier,  1997;  Ramp  et  al.,  2005).  This  lead  to  a  very  strong  gradient 
in  the  wind  stress  extending  southwest  from  Santa  Cruz. 

To  quantify  these  gradients,  the  wind  stress  curl  (V  x  t)  was  computed  on  the 
same  grid  as  the  wind  stress  itself  (Figure  4).  A  planer  least-squares  fit  was  used  over  a 
circle  with  radius  15  km  to  estimate  the  curl.  During  the  weak  poleward  winds  on 
August  4-5,  the  curl  changed  sign  down  the  middle  of  the  observing  region  about  20  km 
offshore  and  was  positive  (upwelling  favorable)  offshore  and  negative  (downwelling 
favorable)  nearshore.  This  may  have  contributed  to  the  rapid  wanning  of  the  sea  surface 
(>3°C)  between  the  4th  and  the  5th  off  Ano  Nuevo.  On  August  8,  10,  and  11,  there  was  a 
region  of  high  positive  curl  extending  southward  across  the  Monterey  Bay  from  Santa 
Cruz,  along  the  boundary  between  the  unobstructed  offshore  wind  and  the  wind  shadow 
in  the  inner  bay.  Using  the  conventional  fonnula  for  F.kman  pumping 

we  =  (V  x  x)  /  pf 
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the  upwelling  velocity  we  was  computed  where  p  is  the  water  density,  f  is  the  Coriolis 
parameter,  and  the  observed  curl  maximum  =  6x10'  N  m'  .  The  result  =  5.8  m  d'  and 
was  comparable  to  the  upwelling  due  to  coastal  divergence  at  the  same  time,  which  was 
in  the  range  of  5-10  m  d’1.  This  indicates  that  the  cold  water  extending  southward  across 
the  Monterey  Bay  during  upwelling  events  may  result  from  the  local  wind  stress  curl 
maximum  in  addition  to  the  southward  advection  of  cold  water  from  the  north.  This  is 
consistent  with  previous  results  on  larger  spatial  scales  that  also  show  that  the  vertical 
velocities  due  to  Ekman  pumping  are  comparable  to  those  from  Ekman  transport  near 
headlands  along  the  central  California  coast  [Pickett  and  Paduan,  2003]. 

3.2  ADCP  Currents 

The  most  obvious  feature  in  the  observed  currents  was  the  preponderance  of 
poleward  flow  at  both  locations  (Figure  2).  One  particularly  strong  poleward  event 
started  on  August  3rd  at  ADCP  2  and  August  4  at  ADCP  1.  The  event  was  nearly 
barotropic  at  both  moorings  with  a  slight  mid-depth  maximum  (near  40  cm  s'1),  and 
ended  around  August  7.  The  temporal  agreement  with  the  August  4-7  wind  relaxation 
appears  coincidental,  in  that  poleward  flows  with  this  magnitude  did  not  appear  during 
the  subsequent  wind  relaxation  events.  Plan-view  displays  of  the  model  output,  discussed 
in  detail  in  a  later  section,  gives  some  indication  that  this  strong  poleward  event  was 
confined  to  the  continental  shelf  and  was  associated  with  a  pulse  of  warm  water  exiting 
the  Monterey  Bay.  There  was  likewise  evidence  in  the  models  that  the  equatorward  flow 
at  ADCP  2  centered  on  August  11-12  was  due  to  this  instrument  being  located  in  the 
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coastal  upwelling  jet,  while  ADCP  1  was  located  well  shoreward  of  the  jet  where  the 
flow  remained  poleward. 

Two  subsequent,  weaker  poleward  flow  events  took  place  from  August  19-21  and 
August  26-28  with  no  obvious  phase  lag  between  moorings.  These  events  were 
barotropic  in  nature  but  had  stronger  currents  near  the  bottom  (order  15  cm  s'1)  than  the 
surface  (order  5-10  cm  s'1).  These  poleward  flows  took  place  during  periods  of  no  winds 
or  weak  equatorward  winds  and  thus  appear  to  be  wind-related.  The  vertical  current 
profdes  during  these  times  could  be  modeled  as  the  sum  of  a  barotropic  poleward 
pressure  gradient  force  over  the  entire  water  column  and  weak  equatorward  Ekman 
forcing  in  the  upper  20  m.  This  Ekman  flow  was  clearly  observed  in  the  across-shore 
component  as  offshore  flow  (blue  colors  in  Figure  2)  during  the  stronger  upwelling- 
favorable  wind  events,  as  for  instance  during  August  10-11,  16-17,  22-25,  and  29-31, 
most  clearly  at  ADCP  2.  In  general,  ADCP  2  showed  more  visual  correlation  with  the 
surface  wind  stress  than  ADCP  1,  and  this  was  supported  by  the  spectral  analysis  (not 
shown).  Significant  spectral  peaks  were  observed  for  periods  of  10  days,  5  days,  and  the 
tides.  A  coherence  and  phase  analysis  showed  that,  on  average,  the  5-10  d  winds  lead  the 
alongshore  currents  by  17.9  hours  at  both  instruments.  This  is  close  to  the  local  inertial 
period  (20  hours)  and  represents  the  spin-up  time  for  the  alongshore  flow  under  wind 
forcing.  At  ADCP  2,  the  across-shore  currents  at  the  near-surface  and  near-bottom  were 
coherent  with  the  alongshore  wind  stress  with  opposite  phase  (151°),  consistent  with 
Ekman  dynamics.  At  ADCP  1,  there  was  no  coherence  between  alongshore  winds  and 
across-shore  currents  at  any  depth.  This  may  be  due  to  the  coastal  constraint,  as  Ekman 
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turning  is  reduced  the  closer  the  location  is  to  the  coast  [Ekman,  1905;  Ramp  and  Abbott, 
1998], 

Also  obvious  in  the  plots  are  the  much  shorter  space  and  time  scales  in  the  across- 
shore  flow  vs.  the  along-shore  flow.  The  across-shore  flow  was  dominated  by  features 
order  20  m  vertically  and  about  2  days  in  length.  There  is  some  suggestion,  especially  at 
ADCP  1,  that  these  features  were  propagating  upward.  Alongshore  currents  were 
statistically  well  correlated  between  moorings  (r  =  .72-. 86)  but  across-shore  currents  were 
weakly  correlated  (r  =  .3-. 5).  The  origin  of  these  smaller-scale  cross-shelf  motions  is  not 
obvious  from  this  data  set.  The  smaller-scale  motions  may  have  an  alongshore 
component  as  well,  but  it  is  masked  by  the  larger-scale  alongshore  flows. 


4.  Model  /  Data  Comparisons 

4.1  Basic  description  of  the  three  numerical  models 

A  key  goal  of  the  ASAP  program  was  to  collect  data  sets  to  improved  numerical 
models  for  predicting  ocean  conditions  at  the  scales  necessary  to  resolve  coastal 
upwelling  and  subsequent  plume  dynamics.  The  models  were  run  in  near-real  time  to 
provide  knowledge  input  to  the  adaptive  sampling  decision  tree.  The  three  models  used 
were  the  Harvard  Ocean  Prediction  System  (HOPS)  [Robinson,  1999;  Haley  et  ah,  2009; 
Haley  and  Lermusiaux,  2010],  the  Jet  Propulsion  Laboratory  implementation  of  the 
Regional  Ocean  Modeling  System  (JPL/ROMS)  [Schepetkin  and  McWilliams,  2004], 
and  the  Navy  Coastal  Ocean  Model  /  Innovative  Coastal-Ocean  Observing  Network 
(NCOM/ICON)  model  [Shulman  et  ah,  2002;  2004],  Two  versions  of  the  HOPS  model 
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were  examined:  HOPS  solo  (hereafter  HOPSs)  which  had  1.5  km  horizontal  resolution, 
and  HOPS  nested  (hereafter  HOPSn)  which  was  nested  within  HOPSs  and  had  0.5  km 
horizontal  resolution  (Figure  5).  JPL/ROMS  and  NCOM/ICON  had  rectangular  grid 
cells  near  the  moorings  with  about  1.5  x  2.5  km  horizontal  resolution.  All  the  models  had 
30-32  sigma  coordinates  in  the  vertical.  While  conceptually  similar  in  operation,  the 
models  differed  significantly  in  tenns  of  the  mixing  schemes  employed,  initial  conditions 
used,  and  the  larger-scale  models  they  were  nested  within.  All  four  models  were  forced 
using  surface  wind  stresses  and  heat  fluxes  from  the  Coupled  Ocean/Atmosphere 
Mesoscale  Prediction  System  (COAMPS™),  and  assimilated  the  NPS  aircraft  SST,  glider 
T  and  S,  and  AUV  T  and  S  data.  Additionally,  the  JPL/ROMS  model  assimilated  T  and 
S  data  from  the  MBARI  buoys  Ml  and  M2.  None  of  the  models  assimilated  velocity 
from  anything  including  vehicles,  moorings,  and  coastal  HF  radar.  The  models  all  use 
different  mixing  parameterizations  (see  references  for  details).  The  effect  of  this  on 
mesoscale  model  output  differences  is  not  obvious,  but  is  expected  to  be  limited  on  a  few 
days  of  ocean  forecasts.  Two  factors  that  clearly  are  important  are  the  model  initial  and 
open  boundary  conditions.  Given  the  small  scale  of  the  domain  off  Point  Ano  Nuevo,  the 
boundary  forcing  is  extremely  important  to  obtaining  the  correct  answer  within  the  model 
domain.  Such  well-known  features  as  the  California  Undercurrent  over  the  continental 
slope  and  the  barotropic  poleward  flow  over  the  continental  shelf  were  clearly  remotely- 
forced  and  entered  the  domain  through  the  open  boundaries.  The  JPL/ROMS  model  had 
the  longest  spin-up  time  using  climatology  since  July  2003,  using  boundary  conditions 
from  the  next-larger  ROMS  model  (5-km)  it  was  nested  within.  NCOM/ICON  likewise 
used  boundary  conditions  from  the  next  larger  nest,  called  NCOM  California  Current 
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System  (NCOM/CCS).  The  run  was  initialized  using  more  recent  data  starting  July  11, 
2006.  HOPS  solo  was  unique  among  the  three  in  that  is  was  not  nested  into  a  larger-scale 
model.  This  model  was  initialized  using  local  data,  that  is,  the  early  glider  runs  and  ship 
surveys  from  the  POINT  SUR  and  the  WECOMA  conducted  in  the  area,  and  employed 
modified  Orlanski  radiation  conditions  at  the  boundaries.  All  the  regional  models,  either 
directly  or  indirectly,  made  use  of  satellite  products  such  as  sea  surface  temperature 
(SST)  and  surface  height  (SSH).  Some  of  the  runs  (such  as  5  km  JPL/ROMS) 
assimilated  satellite  data  directly,  while  the  HOPS  runs  did  not  assimilate  SSH. 
NCOM/CCS  ingested  SSH  via  an  intennediate  step,  the  Navy’s  Modular  Ocean  Data 
Assimilation  (MODAS)  system.  This  infonnation  was  thus  transferred  to  the  smaller  1 .0- 
1.5  km  nests  via  the  open  boundary  conditions. 

4.2  Visual  comparisons 

To  develop  a  basis  for  comparison,  the  output  from  all  four  numerical  models  first 
underwent  the  same  low-pass  filtering  and  coordinate  rotation  as  the  ADCP  data.  The 
principal  axis  ellipses  are  interesting  in  their  own  right:  At  ADCP  1  the  NCOM  model 
overestimated  the  currents  by  about  2  cm/s  relative  to  both  the  remaining  models  and  the 
data,  while  the  ROMS  model  output  was  misaligned  by  12  degrees  relative  to  the  other 
estimators  (Figure  6a).  At  ADCP  2,  all  the  models  statistically  underestimated  the  true 
currents  and  the  HOPSs  model  had  too  much  energy  in  the  across-shelf  component 
(Figure  6b). 

Three  depths  were  chosen  for  visual  and  statistical  comparison:  Near-surface  (10, 
12  m),  mid-depth  (24,  52  m),  and  near-bottom  (46,  72  m)  for  ADCPs  (1,2)  respectively. 
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Recall  that  at  ADCP  1,  the  observed  currents  consisted  of  strong  poleward  flow  from  3-8 
August,  followed  by  equatorward  flow  from  9-17  August,  and  weakly  poleward  flow  for 
the  remainder  of  the  record  (18-31  August).  This  was  true  at  all  depths  with  limited 
vertical  shear  (Figure  7).  The  NCOM  model  fluctuations  appeared  to  lead  the  data,  most 
noticeably  near  the  surface  (Figure  7a).  After  mid-August,  the  NCOM  model 
overestimated  the  strength  of  the  poleward  flow  by  a  factor  of  2-4,  worse  near  the  surface 
than  near  the  bottom.  The  HOPSs  model  overestimated  the  poleward  flow  from  1-3 
August,  missed  the  equatorward  event  completely,  and  also  overestimated  the  second 
poleward  event  (Figure  7a).  Like  the  NCOM  model,  the  model/data  agreement  improved 
with  depth  and  near  the  bottom  the  agreement  was  much  better  than  higher  in  the  water 
column  (Figure  7c).  The  ROMS  model  showed  better  agreement  with  the  data  near  the 
surface  than  the  NCOM  or  HOPSs  models  but  still  had  difficulty  reproducing  the 
equatorward  flow.  The  ROMS  currents  opposed  the  observed  currents  during  9-13 
August  and  were  too  weak  from  13-18  August.  The  ROMS  output  was  the  least  sheared 
among  the  models  and  looked  much  the  same  at  all  depths  (Figure  7,  bottom  panels). 
The  HOPSn  model  showed  the  best  visual  agreement  with  the  data.  Like  HOPSs,  the 
HOPSn  time  series  began  with  too  much  poleward  flow,  perhaps  a  product  of  model 
initialization  and  the  short  spin-up  time  relative  to  the  other  two  models.  Once  it  got  on 
track,  the  model  produced  the  August  13-17  equatorward  event  and  had  close  to  the  right 
magnitude  for  the  second  poleward  event. 

The  observations  at  ADCP  2  were  more  vertically  sheared  than  at  ADCP  1 .  At 
ADCP  2,  the  near-surface  currents  were  entirely  equatorward  following  the  wind  stress 
during  August  8-31,  while  the  mid-depth  and  near  bottom  currents  were  poleward  from 
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August  18-28  in  opposition  to  the  surface  flow  (Figure  8,  compare  a  and  b).  The 
observed  surface  currents  at  ADCP  2  also  opposed  the  surface  currents  at  ADCP  1  from 
August  18-28,  indicating  the  presence  of  a  velocity  front  between  the  two  instruments 
(compare  Figures  7a  and  8a).  A  model’s  perfonnance  will  therefore  be  quite  sensitive  to 
its  ability  to  correctly  position  that  front.  Like  at  ADCP  1,  the  NCOM  model  looked 
good  for  the  first  half  of  the  record,  again  leading  the  data  slightly,  but  deteriorated 
during  the  second  half  of  the  month  where  the  surface  poleward  flow  in  the  model  output 
opposed  the  observed  flow  at  the  ADCP  (Figure  8a).  This  was  perhaps  because  aircraft 
and  glider  data  were  available  to  assimilate  during  August  1-15  but  not  during  August  16- 
3 1 .  The  NCOM  model  and  data  agreement  at  ADCP  2  was  much  better  near  the  bottom 
where  the  model  captured  the  poleward/equatorward/poleward  transitions  although  the 
timing  was  still  off,  and  there  was  too  little  equatorward  flow  and  too  much  poleward 
flow  at  the  bottom.  The  equatorward  event  was  weak  in  the  ROMS  model  near  the 
surface  and  absent  at  the  other  depths,  indicating  it  misplaced  the  location  of  the 
California  Current  boundary  (Figure  8a-c).  It  also  badly  overestimated  the  poleward  flow 
after  August  18.  Like  the  NCOM  model,  the  currents  after  August  18  more  closely 
resembled  the  observations  at  ADCP  1  than  at  ADCP  2,  indicating  that  the  model  placed 
the  velocity  front  offshore  of  the  observed  velocity  front.  This  is  consistent  with  both 
NCOM  and  ROMS  currents  being  geographically  misplaced  a  bit  too  far  offshore  of  the 
real  ocean.  The  HOPSs  model  looked  good  at  the  beginning,  but  then  had  equatorward 
rather  than  poleward  flow  at  mid-depth  and  near-bottom  after  August  18  (Figure  8b,  c). 
The  HOPSn  model  looked  the  best  at  site  2  as  it  did  at  site  1 .  It  faithfully  produced  both 
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the  opening  poleward  and  equatorward  flow  events,  and  at  least  had  the  sign  correct  at 
the  end. 


4.3  Statistical  comparisons 

While  already  interesting,  these  stick  vector  plots  form  merely  an  introduction  to 
more  quantitative  statistical  comparisons.  The  tools  used  to  facilitate  these  comparisons 
were  the  temporal  cross-correlation  functions  and  the  root  mean  square  (RMS)  error 
between  the  model  output  and  the  ADCP  data  at  several  grid  points  surrounding  each 
ADCP  location  (Figure  9).  The  correlations  determine  if  the  model  output  was 
geographically  or  temporally  displaced  from  the  data,  while  the  RMS  error  estimates  the 
magnitude  of  the  modcl/data  differences. 


For  the  closest  model  grid  point  to  the  ADCP  location  and  each  of  the 
surrounding  grid  points,  the  correlation  coefficient  p  (z)  was  calculated  as  [Bendat  and 

xy 

Piersol,  1986]: 


where  Cxy  (x)  is  the  cross-covariance  between  the  model  output  (x)  and  the  data  (y), 
crv,  <jv  are  the  variances  of  series  x  and  y,  and  x  is  the  time  lag.  Significance  levels  for 

the  correlations  were  calculated  following  the  method  of  Amos  and  Koopmans  [1963] 
using  the  effective  degrees  of  freedom  estimates,  Nejf,  defined  as  [Davis,  1976]: 

Neff=N/Tn 

where  N  is  the  length  of  the  time  series  in  hours  and 

T.  =ZAC»('A,)C.('A')A< 
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is  the  correlation  time  scale  for  the  data.  More  self-correlated  data  results  in  fewer  true 
degrees  of  freedom  and  therefore  higher  significance  levels  for  the  correlation  coefficient. 

When  examining  the  results  from  these  calculations  it  is  important  to  keep  the 
spatial  distribution  of  the  points  being  compared  in  mind  as  they  were  different  for  each 
model  (Figure  9).  For  the  NCOM  and  both  HOPS  models,  the  grid  points  were  oriented 
along/across-shore  while  the  ROMS  model  grid  was  oriented  north/south.  The  HOPSn 
grid  was  about  3  times  more  resolved  than  the  others  and  less  variation  is  thus  expected 
in  the  correlation  coefficients  between  these  points  and  the  ADCP  at  the  center.  The 
results  (Appendix  A)  are  displayed  as  box  charts  showing  the  correlation,  RMS  error,  and 
optimum  lag  time  for  each  of  the  nine  grid  points  compared  to  the  ADCP  located  at  the 
center.  Results  which  were  statistically  insignificant  are  shaded,  and  time  lags  are 
presented  only  for  the  statistically  significant  results.  A  negative  (positive)  phase  means 
the  model  leads  (lags)  the  data.  An  empty  box  indicates  no  data,  usually  because  the 
selected  model  depth  was  deeper  than  the  bottom  on  the  inshore  side  of  the  mooring.  For 
display,  the  along-  and  across-shore  components  for  each  model  from  the  grid  point 
closest  to  the  ADCP  were  compared  with  the  data  (Figure  10).  Then,  the  tables  were 
used  to  produce  component  plots  of  the  “best”  modePdata  comparisons  for  each  model 
using  the  grid  point  with  the  highest  correlation  shifted  relative  to  the  data  by  the 
corresponding  time  lag  (Figure  1 1).  A  few  robust  results  emerged  from  these 
comparisons. 

For  the  alongshore  component,  all  the  models  except  NCOM  at  ADCP  1  were 
significantly  correlated  with  the  data.  At  site  1,  the  ROMS  model  produced  the  highest 
correlations  followed  by  HOPSs  and  HOPSn.  The  ROMS  model  output  was  very  nearly 
in  phase  with  the  observations  while  both  the  HOPS  models  lead  the  data  by  about  16  - 
24  hours  (Figure  10,  left  column).  At  site  2,  the  highest  correlations  were  at  HOPSn 
followed  closely  by  HOPSs  and  ROMS  where  the  correlation  levels  were  quite  similar. 
The  HOPSn  model  lead  the  data  by  about  5  hours  at  the  surface  but  lagged  by  16  to  21 
hours  subsurface.  The  HOPSs  model  also  lead  the  data  but  by  smaller  amounts,  0-8 
hours  for  the  most  correlated  points.  The  phase  numbers  for  ROMS  at  site  1  were  a  bit 
erratic:  Choosing  the  most  correlated  points,  the  model  was  in  phase  as  the  surface,  lead 
the  data  by  22  hours  at  mid-depth,  and  lagged  by  7  hours  at  the  bottom  (Appendix  A).  In 


21 


terms  of  the  geographic  accuracy,  HOPSn  at  site  2  was  consistently  better  correlated  at 
the  NE  point  rather  than  at  the  center  point,  perhaps  indicating  a  slight  misalignment  of 
model  features,  although  the  distances  were  small.  None  of  the  other  models  did 
consistently  better  or  worse  at  the  other  model  grid  points. 

The  HOPS  nested  model  was  the  only  one  with  any  skill  at  reproducing  the 
across-shelf  component.  Some  significant  correlations  were  observed  for  the  ROMS 
model  but  they  were  180  out  of  phase.  At  site  2,  the  correlations  between  HOPSn  and  the 
data  were  well  above  the  significance  levels  (Appendix  A)  and  the  visual  agreement  with 
the  observations  was  good  (not  shown).  This  indicates  that  HOPSn  responded  well  to  the 
local  wind  forcing  as  described  earlier  from  the  observations  alone,  which  was  more 
evident  at  site  2  than  at  site  1 .  The  HOPSn  across-shore  component  always  lagged 
behind  the  observations  by  longer  (18  hours)  at  mid-depth  than  at  the  surface  and  bottom 
(8  hours).  This  is  again  consistent  with  a  boundary  layer  response,  however  it  takes  too 
long  in  the  model  for  the  momentum  to  be  transferred  to  mid-depth.  It  seems  that  very 
high  resolution  models  are  necessary  to  successfully  simulate  the  across-shore 
component  of  the  flow  on  the  continental  shelf  in  this  region.  This  is  likely  due  to  the 
shorter  correlation  scales  in  the  across-shelf  direction  than  the  alongshelf  direction  [Ramp 
et  ah,  2008]. 


5.  Discussion 

5.1  Remote  Forcing  and  Boundary  Conditions 

All  three  models,  with  well-established  performance  at  larger  space  and  time 
scales,  had  difficulty  reproducing  the  current  variability  in  this  small  sample  region  off 
Point  Ano  Nuevo.  One  potential  explanation  is  that  these  small  domains  cannot  correctly 
capture  remote  forcing  in  the  form  of  alongshore  pressure  gradient  forces  or  coastally- 
trapped  waves,  which  propagate  from  south  to  north  with  the  coast  on  the  right  in  this 
region.  This  means  that  larger-scale  features  such  as  these  must  be  fed  to  the  smaller 
regional-scale  models  via  their  boundary  conditions.  Two  of  the  models,  ROMS  and 
NCOM,  obtained  these  boundary  conditions  from  the  models  they  were  nested  within. 
The  boundaries  in  HOPSn  were  obtained  from  HOPSs,  which  was  initialized  using  the 
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ASAP  opening  ship  surveys  and  other  available  data.  The  HOPSs  boundary  conditions 
were  subsequently  updated  using  fits  of  the  model  to  all  available  data. 

To  explore  the  effects  of  different  boundary  conditions  and  coordinate  systems, 
additional  model  runs  were  made  subsequent  to  the  ASAP  program  using  the  Navy’s 
global  Hybrid  Coordinate  Ocean  Model  (HYCOM)  to  force  the  NCOM/ICON  model, 
rather  than  the  global  NCOM  via  the  NCOM  CCS  output.  Forcing  the  same  NCOM 
ICON  model  directly  from  the  global  HYCOM  model  produced  a  dramatic  improvement 
in  the  model/data  comparisons:  The  complex  correlations/angular  displacements  (Kundu, 
1976)  improved  from  0.04/52.0°  to  0.69/-1 1.5°  at  ADCP  1,  and  from  0.17/2.9°  to  0.86/- 
4.9°  at  ADCP  2.  This  difference  represents  a  change  from  essentially  no  predictive  skill 
to  very  good  predictive  skill. 

These  improvements  can  be  understood  by  comparing  the  sea  surface  height 
(SSH)  anomalies  along  the  California  coast  calculated  from  the  NCOM  /  NCOM  CCS 
model  and  the  HYCOM  model  with  the  coastal  SSH  observations  (Figure  12).  Both 
models  and  the  observations  show  that  the  SSH  anomalies  during  August  2006  were 
larger  than  normal  and  much  larger  than  during  AOSN-II  in  summer  2003  [Ramp  et  ah, 
2009].  Both  models  also  show  that  these  SSH  anomalies  were  propagating  poleward  at 
1-3  m  s'1,  as  can  bee  seen  by  the  slope  of  the  maxima  (red  region)  in  the  latitude  vs.  time 
plots  (Figure  12).  These  anomalies  therefore  represent  remote  forcing  which  will  impact 
the  smaller  ASAP  region  off  Point  Ano  Nuevo.  Comparing  the  model  SSH  (red,  blue 
lines)  with  the  observed  SSH  at  Monterey  and  San  Diego  (black  lines)  shows  that  NCOM 
/  NCOM  CCS  model  underestimated  these  propagating  anomalies  with  respect  to  both 
the  observations  and  the  HYCOM  model,  which  closely  tracked  the  observed 
fluctuations. 

Many  previous  authors  have  documented  poleward-propagating  motions  along  the 
U.S.  west  coast  in  coastal  sea  level  records  [e.g.  Chelton  and  Davis,  1982;  Denbo  and 
Allen,  1987;  Spillane  et  al.,  1987]  current  meter  records  [Ramp  et  ah,  1997a]  and  model 
output  [Pares-Sierra  and  O’Brien,  1989;  Ramp  et  al.  1997b].  For  ASAP  2006,  a 
coherence  and  phase  analysis  of  the  observed  SSH  anomalies  relative  to  the  annual  mean 
along  the  coast  during  May  to  October  2006  showed  poleward  propagation  at  2.7  m  s'1 
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between  Port  San  Luis  and  Monterey  for  oscillations  with  a  15-day  period.  Using  a 
simple  c  =  co/k  relationship  where  c  is  the  phase  speed,  co  is  the  angular  frequency,  and  k 
is  the  horizontal  wave  number,  these  fluctuations  would  have  a  wavelength  of  3695  km  if 
they  were  propagating  as  free  waves  along  the  coast.  Previous  authors  found  waves 
along  the  Big  Sur  coast  in  the  29-day  band  that  closely  resembled  linear,  coastally- 
trapped  mode-1  internal  Kelvin  waves  [Ramp  et  ah,  1997a],  The  dispersion  relation  for 
these  waves  is  given  by 

2  _  N2k2 
(°n  ~  k 2  +(nj t/H)2 

where  N  is  the  Brunt-Vaisala  frequency,  n  is  the  mode  number,  H  is  the  bottom  depth, 
and  the  other  variables  are  as  previously  defined.  Choosing  N  =  5.5  x  10'  s'  as  a 
representative  value  and  H  =  1500  m  as  a  mean  depth  over  the  slope,  the  dispersion 
relation  gives  a  wavelength  of  3490  km  for  the  15-day  period.  The  fluctuations  observed 
in  both  coastal  sea  level  and  the  Navy  models  are  therefore  consistent  with  coastally- 
trapped  mode-1  internal  Kelvin  waves.  It  appears  likely  therefore  that  the  improvement 
in  the  model/data  comparisons  at  ADCP  1  and  ADCP  2  with  HYCOM  nesting  were  due 
to  that  model’s  more  accurate  representation  of  the  poleward  fluctuations  entering  the 
ASAP  region  from  the  south.  Additional  model  runs  (not  shown)  changing  only 
smoothed  (sigma)  vs.  un-smoothed  (hybrid)  bottom  coordinates  showed  improvements 
using  the  hybrid  coordinates,  but  these  improvements  were  much  smaller  than  those 
resulting  from  using  the  HYCOM  vs.  NCOM  model.  This  indicates  that  the  improved 
boundary  conditions  were  the  most  important  factor  in  improving  the  NCOM  ICON 
model  perfonnance. 


5.2  Importance  of  Offshore  Eddies 

Given  the  relatively  good  visual  and  statistical  comparisons,  a  few  plan-view  plots 
from  the  HOPSn  model  were  chosen  to  illuminate  the  larger-scale  context  for  the  ADCP 
observed  currents  (Figure  13).  The  first  pair  on  August  4  and  6  (Figures  13  a  and  b) 
show  the  development  of  the  strong  poleward  flow,  one  of  the  most  prominent  features  of 
the  record.  The  cold  water  at  Point  Ano  Nuevo  was  quickly  swamped  by  the  wann  water 
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from  the  south  during  this  event.  The  only  really  strong  equatorward  flow  during  August 
2006  took  place  during  the  days  surrounding  August  1 1 .  This  was  due  to  the 
development  of  a  coastal  upwelling  jet  that  was  sandwiched  between  the  upwelling  center 
and  a  warm  shoreward  incursion  of  the  California  Current  offshore  (Figure  13c).  During 
subsequent  upwelling  events,  as  exemplified  for  instance  by  conditions  on  August  23 
(Figure  13d),  the  upwelling  front  was  located  further  offshore  and  the  equatorward  flow 
at  the  two  ADCPs  was  weaker  or  non-existent.  A  comparison  of  Figures  13a  and  13c 
supports  the  conclusion  of  earlier  investigators  that  the  bifurcation  of  the  upwelling 
plume  is  closely  related  to  the  position  and  strength  of  the  offshore  eddy  field.  When  the 
eddy  was  well  offshore  (Figure  13a)  the  cold  plume  likewise  advected  offshore.  When 
the  eddy  was  strong  and  closer  to  shore  (Figure  13c)  the  cold  plume  moved  primarily 
southward. 


6,  Summary  and  Conclusions 

The  ocean’s  response  to  both  upwelling  and  relaxation  events  was  studied  for  the 
month  of  August  2006  near  Point  Ano  Nuevo,  CA  as  part  of  the  ASAP  MB06  field 
experiment.  An  aircraft,  several  autonomous  gliders  and  AUVS,  moored  ADCPs,  and 
three  numerical  models  were  used  to  describe  the  flow.  This  paper  provides  a  basic 
description  of  the  flow  from  the  aircraft  and  moorings  and  compares  those  results  with 
numerical  simulations  from  the  HOPS,  JPL/ROMS,  and  NCOM  models. 

The  aircraft  produced  a  high-resolution  spatial  description  of  the  wind  stress 
variability  which  was  previously  unavailable  for  the  region.  Several  new  features 
emerged  including  the  wind  stress  maximum  to  the  north  of  the  Point  during  poleward 
wind  events,  a  clear  maximum  offshore  from  the  Point  during  upwelling  events,  and  a 
strong  wind  stress  curl,  also  during  upwelling  favorable  events,  which  could  potentially 
drive  upwelling  with  strength  (5-10  m  d'1)  equal  to  the  coastal  divergence.  This  spatial 
variability  in  the  wind  itself  likely  plays  an  important  role  in  determining  the  location  of 
3-D  upwelling  centers  along  the  west  coast  of  the  United  States.  The  wind  stress 
observed  at  the  closest  local  buoy  was  often  not  representative  of  the  actual  wind  stress 
off  the  Point  or  to  the  north  of  the  Point  and  much  caution  is  needed  when  applying  the 
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buoy  wind  observations  to  the  local  coastal  oceanography.  The  MBARI  buoy  M2  was 
nevertheless  the  most  representative  time  series  available  for  the  entire  month  of  August 
since  the  aircraft  was  not  available  after  August  15.  The  buoy  recorded  upwelling- 
favorable  wind  events  with  approximately  equal  strength  during  August  1-3,  8-12,  16-25, 
and  28-31.  Poleward  winds  less  than  2.5  m  s'1  blew  during  August  4-7,  13-15,  and  26- 
27. 

Despite  being  only  6.5  km  apart,  moorings  1  and  2  responded  differently  to  the 
sequence  of  upwelling  favorable  winds  separated  by  brief  relaxations.  The  flow  at  site  1 
was  almost  always  poleward  and  often  appeared  to  surge  out  of  the  Monterey  Bay.  Site  2 
was  just  far  enough  offshore  to  miss  this  phenomena.  The  flow  there  was  more  strongly 
sheared  vertically  and  more  coherent  with  the  local  wind  stress,  especially  the  across- 
shore  current  component  with  the  alongshore  wind.  Still,  there  were  many  times  when 
the  currents  at  site  2  were  not  coherent  with  the  winds  and  were  instead  well  correlated 
with  the  position  of  the  mesoscale  features  in  the  coastal  ocean.  The  strongest 
equatorward  currents  were  clearly  associated  with  the  upwelling  jet,  and  the  across-shore 
position  of  the  jet  largely  determined  the  current  strength  at  site  2,  less  so  at  site  1  which 
was  located  well  inshore  of  the  jet  most  of  the  time. 

Given  the  very  small  scales  of  motion  and  the  strong  shears  between  the  two 
moorings,  it  is  perhaps  not  surprising  that  the  three  models  all  had  difficulty  simulating 
the  flow.  The  models  had  more  success  during  the  first  half  of  the  month  when  the  flow 
was  more  wind-driven  and  there  were  aircraft  data  available  to  assimilate.  The  models 
struggled  during  the  second  half  of  the  experiment  when  the  flow  appeared  to  be  more 
remotely-forced  and  local  observations  were  lacking.  A  post-experiment  analysis  of 
observed  sea  level  anomalies  along  the  coast  and  output  from  the  global  HYCOM  model 
re-analysis  both  showed  that  sea  level  fluctuations  during  August  2006  were  anomalously 
strong  and  propagated  poleward  with  characteristics  consistent  with  free  coastally- 
trapped  mode-1  internal  Kelvin  waves.  The  failure  of  the  outer-nests  driving  the  open 
boundaries  of  the  JPL-ROMS  and  NCOM  ICON  models  to  correctly  capture  these 
fluctuations  may  account  for  many  of  the  observed  weaknesses  in  the  model/data 
comparisons.  Subsequent  re-analyses  of  the  NCOM  ICON  model  forced  by  the  global 
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HYCOM  model,  which  more  accurately  represented  the  coastal  sea  level  variability, 
showed  a  dramatic  improvement  in  model  predictive  skill. 

The  only  model  that  had  any  demonstrated  skill  simulating  the  across-shore 
component  was  the  HOPS  nested  model  (HOPSn),  which  had  a  three  times  finer  grid 
resolution  than  the  others  (0.5  vs.  1.5  km).  This  very  fine  simulation  scale  appears  to  be 
necessary  to  correctly  place  the  strong  fronts  that  often  occur  in  the  coastal  ocean. 
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Figure  2.  Time  series  of  a)  wind  stress  from  MBARI  Buoy  M2,  b)  across-shore  current 
at  ADCP  1,  c)  alongshore  current  at  ADCP  1,  d)  across-shore  current  at  ADCP  2,  and  e) 
alongshore  current  at  ADCP  2.  The  scale  of  the  y-axis  is  the  same  for  both  ADCPs,  with 
the  difference  in  plot  height  reflecting  the  difference  in  water  depth.  The  velocity  scales 
are  different  for  the  across-  and  alongshore  components  to  better  illustrate  the  much 
smaller  across-shore  flow.  A  positive  sign  indicates  poleward  and  onshore  flow.  The 
vertical  gray  bars  on  the  wind  stress  plot  indicate  the  time  of  the  TWIN  OTTER  over¬ 
flights. 
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Figure  3.  Time  series  of  the  aircraft-observed  surface  wind  stress  during  (a-d)  the  August 
2-5  wind  relaxation  event  and  (e-h)  the  August  8-12  upwelling  favorable  event. 
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Figure  4.  The  curl  of  the  wind  stress  computed  on  the  same  grid  as  Figure  3.  The  units 

■5 

on  the  color  bar  are  in  N  m'  . 
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Longitude  (deg  W) 

Figure  5.  Model  domains  off  central  California  for  the  four  numerical  models  compared 
against  the  moored  ADCP  data  on  the  continental  shelf.  The  mooring  positions  are 
indicated  by  diamonds. 
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Figure  6.  Mid-water  column  principal  axis  ellipses  for  both  models  and  moorings  from: 
A)  ADCP  1,34m  depth;  and  B)  ADCP  2,  54  m  depth.  North  is  towards  the  top  of  the 
plot. 
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Figure  7.  A  comparison  of  currents  from  ADCP  1  and  the  four  numerical  models  for:  a) 
near-surface  (10  m),  b)  mid-depth  (24  m),  and  c)  near-bottom  (46  m).  Poleward 
alongshore  is  up. 
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Figure  8.  A  comparison  of  currents  from  ADCP  2  and  the  four  numerical  models  for:  a) 
near-surface  (12  m),  b)  mid-depth  (52  m),  and  c)  near-bottom  (72  m). 
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Figure  9.  Grid  point  locations  used  to  facilitate  the  spatial  correlations,  with  the  ADCP 
located  at  the  center  of  each  grid  as  indicated  by  the  asterisk. 
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Figure  10.  Model/data  (gray/black  lines)  comparisons  for  the  alongshore  current 
components  at  ADCP  1  (left  column)  and  ADCP  2  (right  column).  The  three  vertical 
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panels  for  both  columns  are  the  comparisons  for  a)  near-surface  (10,  12  m),  b)  mid-depth 
(24,  52  m),  and  c)  near-bottom  (46,  72  m). 
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Figure  11.  As  in  Figure  10,  except  the  plots  are  for  the  grid  points  (see  Figure  9)  and 
time  lags  that  resulted  in  the  best  agreement  between  the  model  output  and  the  ADCP 


HOPSn 

~  HOPSn 

WLag  =  -18hrs  ‘ 

20 

W  Lag  =  -3  hrs 

data,  as  judged  by  the  highest  correlation  coefficient  (Appendix  A).  The  notation 
indicating  the  position  of  the  winning  grid  point  relative  to  the  ADCP  is  the  same  as 
points  on  a  compass  with  a  “C”  indicating  the  center  point  (no  geographical  offset).  The 
convention  for  the  phase  is  that  a  negative  (positive)  sign  means  the  model  leads  (lags) 
the  data.  For  example  (upper  left  panel),  the  alongshore  component  of  the  ROMS  model 
at  10m  showed  the  best  agreement  with  the  data  when  the  southern  grid  point  was  lagged 
6  hours  in  time.  The  gray  shading  indicates  panels  where  the  correlation  coefficient  was 
not  statistically  significant,  but  the  phase  shift  was  used  anyway  for  completeness. 
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Figure  12.  Sea  surface  height  (SSH)  anomalies  from  24°N  to  40°N  as  computed  by  (left 
panel)  the  Navy  regional  NCOM  CCS  model  with  global  NCOM  boundary  conditions 
and  (right  panel)  the  global  HYCOM  model.  The  observed  SSH  from  coastal  sea  level 
observations  at  Monterey  (36°  36'  N)  and  San  Diego  (32°  43'  N)  are  included  as  the 
heavy  black  lines.  Time  series  at  the  same  locations  sub-sampled  from  the  two  different 
model  configurations  are  shown  as  the  red  (NCOM)  and  blue  (HYCOM)  lines 
respectively. 
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Figure  13.  HQPSn  model  currents  (black  vectors)  overlaid  on  temperature  (color  bar)  for 
the  sea  surface  layer  off  Point  Ano  Nuevo  for  a)  4  August,  b)  6  August,  c)  1 1  August,  and 
d)  23  August  2006. 
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APPENDIX  A. 


Statistical  comparisons  of  the  model  output  at  nine  grid  points  with  the  ADCP  time 
series,  co-located  at  the  center  of  the  box.  The  maximum  correlation  coefficient  at  the 
optimum  lag  (hours),  and  RMS  error  between  instruments  are  shown  for  each 
comparison.  Correlations  that  are  not  statistically  significant  have  been  grayed  out.  The 
across-  and  alongshore  components  are  shown  in  the  left  and  right  column  respectively. 
The  surface,  mid-depth,  and  near-bottom  time  series  are  at  the  top,  middle,  and  bottom. 
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0.28 

Correlation 

0.78 

0.75 

0.73 

RMS  Error 

3.77 

4.88 

4.15 

RMS  Error 

13.75 

7.26 

7.94 

Lag  time  (hrs) 

6 

Lag  time  (hrs) 

-3 

-6 

-15 

0.34 

0.30 

0.20 

0.81 

0.76 

0.73 

8.11 

4.24 

6.91 

15.5224 

8.6336 

7.9259 

-3 

7 

-7 

0.22 

0.21 

0.15 

0.74 

0.76 

0.73 

8.13 

3.93 

7.05 

13.61 

12.44 

9.88 

8 

7 

3 
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HOPS  Nested 
ADCP  1 


Surface  Currents  10.25m 


Across-shelf  Components 


Alongshelf  Components 


Correlation 

0.31 

0.27 

0.21 

Correlation 

0.71 

0.70 

0.69 

RMS  Error 

4.70 

4.52 

4.39 

RMS  Error 

6.44 

6.55 

6.66 

Lag  time  (hrs) 

Lag  time  (hrs) 

-12 

-12 

-12 

0.34 

0.31 

0.26 

0.71 

0.70 

0.69 

4.96 

4.8305 

4.62 

6.34 

6.49 

6.68 

-12 

-12 

-11 

0.36 

0.33 

0.29 

0.71 

0.70 

0.69 

5.04 

5.00 

4.84 

6.18 

6.42 

6.67 

-13 

-12 

-12 

Mid  Depth  34.25m 
Across-shelf  Components 


Alongshelf  Components 


Correlation 

0.55 

0.56 

0.56 

Correlation 

0.69 

0.69 

0.689 

RMS  Error 

1.40 

1.40 

1.51 

RMS  Error 

5.46 

5.47 

5.43 

Lag  time  (hrs) 

10 

8 

8 

Lag  time  (hrs) 

-20 

-21 

-21 

0.53 

0.55 

0.55 

0.69 

0.68 

0.68 

1.49 

1.43 

1.55 

5.45 

5.50 

5.49 

10 

8 

8 

-20 

-20 

-21 

0.51 

0.53 

0.55 

0.69 

0.68 

0.67 

1.60 

1.54 

1.58 

5.42 

5.52 

5.53 

11 

9 

8 

-20 

-20 

-20 

Near  Bottom  46.25m 
Across-shelf  Components 


Alongshelf  Components 


Correlation 

0.26 

0.27 

0.26 

Correlation 

0.66 

0.63 

0.59 

RMS  Error 

1.49 

1.72 

1.92 

RMS  Error 

5.08 

5.27 

5.78 

Lag  time  (hrs) 

Lag  time  (hrs) 

-22 

-21 

-21 

0.24 

0.25 

0.25 

0.66 

0.64 

0.60 

1.52 

1.65 

1.94 

5.07 

5.20 

5.56 

-21 

-21 

-20 

0.25 

0.24 

0.23 

0.67 

0.64 

0.61 

1.57 

1.62 

2.03 

5.02 

5.18 

5.45 

-20 

-21 

-20 
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HOPS  Nested 
ADCP2 


Surface  Currents  12.24m 


Across-shelf  Components 


Alongshelf  Components 


Correlation 

0.64 

0.62 

0.60 

Correlation 

0.88 

0.89 

0.90 

RMS  Error 

5.13 

5.15 

4.88 

RMS  Error 

8.14 

8.36 

8.69 

Lag  time  (hrs) 

7 

8 

9 

Lag  time  (hrs) 

-6 

-5 

-4 

0.62 

0.61 

0.59 

0.87 

0.88 

0.89 

5.14 

5.25 

4.91 

8.08 

8.34 

8.70 

8 

8 

9 

-5 

-5 

-5 

0.61 

0.60 

0.57 

0.86 

0.87 

0.89 

5.14 

5.24 

4.91 

8.09 

8.43 

8.83 

9 

9 

9 

-5 

-5 

-5 

Mid  Depth  52.24m 
Across-shelf  Components 


Alongshelf  Components 


Correlation 

0.52 

0.59 

0.63 

Correlation 

0.85 

0.85 

0.86 

RMS  Error 

1.94 

1.81 

1.68 

RMS  Error 

7.83 

7.77 

7.78 

Lag  time  (hrs) 

17 

17 

17 

Lag  time  (hrs) 

-18 

-19 

-21 

0.52 

0.61 

0.64 

0.84 

0.84 

0.85 

1.92 

1.85 

1.67 

7.84 

7.83 

7.81 

17 

18 

18 

-18 

-19 

-21 

0.50 

0.58 

0.62 

0.84 

0.84 

0.85 

1.93 

1.91 

1.74 

7.78 

7.83 

7.84 

16 

18 

18 

-18 

-19 

-21 

Near  Bottom  72.24m 
Across-shelf  Components 


Alongshelf  Components 


Correlation 

0.45 

0.42 

0.39 

Correlation 

0.81 

0.84 

0.86 

RMS  Error 

1.93 

1.72 

1.85 

RMS  Error 

7.02 

6.67 

6.42 

Lag  time  (hrs) 

11 

3 

Lag  time  (hrs) 

-13 

-16 

-19 

0.48 

0.44 

0.41 

0.81 

0.83 

0.86 

2.15 

1.70 

1.75 

7.25 

6.82 

6.45 

14 

5 

-13 

-16 

-18 

0.51 

0.42 

0.41 

0.80 

0.82 

0.85 

2.20 

1.82 

1.70 

7.43 

7.07 

6.61 

17 

-14 

-16 

-17 
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